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1. Introduction 

Among the few aspects of quantum field theories (QFTs) often amenable to analytical 
study is their behaviour at large and/or short distances. The asymptotic behaviour in 
these regimes is governed by (possibly trivial) fixed points of the renormalization group 
(RG), and RG-improved perturbation theory can provide important insights into the full 
non-perturbative behaviour of the theory. 

A particularly interesting situation arises if both the ultraviolet (UV) and infrared 
(IR) fixed points of a QFT are nontrivial. This is interesting, from a QFT point of view 
because RG flows between such fixed points provide some understanding of the topology of 
the space of QFTs, and from a statistical mechanics viewpoint because it yields examples 
of the cross-over between different universality classes of critical phenomena. One would 
like to understand such RG flows in detail, e.g. not just between which theories the flow 
interpolates but also how specific operators evolve from the UV to the IR. 

In 1+1 dimensions (where the subject may have implications also for the non- 
perturbative formulation of string theory) it is possible to study some of these questions 
explicitly. The best known examples of RG flows between non-trivial fixed points are 
those [0[0 induced by <f>i ^-perturbations of the unitary 0] minimal models J|] of central 
charge c m = 1 — m ^ +1 - ) with diagonal modular invariant partition functions (MIPFs) 0. 
Here we study the c^i^-induced flows between the unitary minimal models with non- 
diagonal MIPFs. 

Denote a unitary minimal model of central charge c m and given MIPF in the ADE 
classification || by X m , with m = 3,4,5,... for X = A, m = 5, 6, 7, . . . for X = D, 
and m = 11,12,17,18,29,30 for X = E. Recall that unlike the A-models, the D- and 
i?-models have non-diagonal MIPFs. Also, except for the .E-models with m = 11, 17, 29 all 
the models contain ^1,3, which is the least relevant (spinless, primary) field.0 Hence, except 
for the latter three models, one can define the perturbed CFTs Xm via the (euclidean) 
action 

A c±) = A x +X I d 2 x 1|3 (a;) , (1.1) 

1 -D5, the 3-state Potts CFT H, is special in that its spectrum contains two copies of ^1,3 
that can be taken to be even (<^ 3 ) and odd (<^ 3 ) with respect to the global Z2 symmetry of the 
model; by "the ^1,3 perturbation" in this case we will always mean the Z2-even perturbation. 
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where the superscripts (±) refer to the sign of A. The theories Xm are integrable 
non-scale-invariant theories, the mass scale being proportional to A 1//ym where y m = 
^q-j- is the RG eigenvalue of the field ^3 (its scaling dimension is d^j = 2 — y m ). 
Since under repeated fusions (^3 closes on itself, modulo irrelevant operators and the 
identity, renomalization of the UV divergences in the perturbation theory based on fll.lQ 
will not generate additional counterterms (the counterterm corresponding to the identity 
does, however, affect the perturbative expansions we will be interested in — see sect. 2). 

Perturbations in opposite directions of the same CFT may lead to completely different 
non-scale-invariant QFTs. A\n , for example, is believed to be a theory of massive kinks ||. 
On the other hand, Am (m > 4) is believed to be massless and flow to A m _\ in the IR, 
as strongly suggested by the LG approach together with perturbative RG analysis 

(the latter applicable when m ^> 1).@ We will present concrete evidence that the theories 
Dm describe the flows 

Dm - ( n 4 'i , (1-2) 

\D m -i if m>6 ' v ; 

the UV (IR) limits corresponding to A = (00), respectively. It is also natural to conjecture 
that Em , Tn = 12, 18, 30, describe the flows E m — > E m -x (see sect. 3.4). 

The two techniques traditionally used to study these and other (conjectured) RG flows 
are Landau-Ginzburg (LG) analysis |nj |T]J [[L2| and RG-improved perturbative calcula- 
tions Hlfl[p||TlT||13|j. Although LG descriptions involving two scalar (real) fields have been 
proposed || Jl4j for D m and E m , they do not seem to be useful when analyzing <pi ^-induced 
RG flows. This is not surprising since the LG models involve two strongly interacting fields, 
even for large m. On the other hand, RG-improved perturbative calculations |[[@ can 
be performed for any provided m ^> 1 (which excludes the ^-models from consid- 

eration), so that the IR fixed point is close to the UV one. This will be discussed in 
sect. 3.3. 

We will also provide evidence for (|1.2|) — in particular for small m — along quite 
different lines. This is possible due to some recent developments in the study of the 
finite-volume spectrum of non-scale-invariant integrable QFTs in 1 + 1 dimensions. The 
thermodynamic Bethe Ansatz (TBA) technique allows one to obtain non-linear integral 



2 The existence of these flows was first alluded to by Huse Q , who noticed a cross-over between 
different critical behaviours in regime IV of the integrable RSOS lattice models underlying the 
QFTs A£\ 
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equations for the exact ground state energy E (R) of such a QFT in finite volume R, 



given its factorizable [15| S-matrix. The small R behaviour of Eq(R) contains information 
about the UV CFT, which can be compared with predictions of conformal perturbation 
theory (CPT). This approach has been extensively used in the last two years as a means 
of checking the purely massive scattering theories conjectured to describe the on-shell 
behaviour of certain perturbed CFTs (see for a brief review of all this, |17] [18| [|19| for 
the details in the case of perturbed CFTs with diagonal S-matrices, and [^D[ for the first 
analysis of a non-diagonal scattering theory). 

Recently [^T| Al. Zamolodchikov proposed an explicit scattering theory of massless 
particles for the non-scale-invariant QFT A^~ , presumably []ll | interpolating between the 
tricritical and the critical Ising CFT. Comparison of the solutions of the TBA equations 
with CPT around the UV and IR CFTs leaves little doubt that the conjectured S'-matrix 
is correct. In [BT| it was further conjectured that Eq(R) in for m > 4 is given by a 



certain generalization (see sect. 2) of the TBA equations for . It was checked that the 
conjectured TBA equations give the expected UV and IR central charges. Nevertheless, 
since for m > 4 the TBA equations have simply been guessed (there are not even any 
conjectures for the S'-matrices of Am^ for m > 4), further evidence for the correctness of 
the equations is required. 

In [B2[ we provided such evidence for m = 5, 6, 7.@ Furthermore, we proposed TBA- 
like equations for certain finite-volume excitation energies in A { ± ] for m even. For A^ 
these equations explicitly demonstrate the flow of the UV conformal fields (f>2,2 and 4>^.,^ 
to fields with the same Kac indices in A m _i, in agreement with expectations based on LG 
and perturbative RG analyses. These results provide further strong (and independent) 
support for the existence of the flows concerned, and show the usefulness of the study of 
finite- volume spectrum in the context of RG flows in general. 

In fact, the "TBA approach" has led the authors of |2c| to conjecture new RG flows 
that has not been proposed before, namely between the Zjv-parafermion CFTs [04] (N > 3) 
and the minimal models A N +\. The conjectures are based on equations for Eq(R) of the 
perturbed parafermion models (the perturbation is by a certain field that breaks the Z N 
symmetry, leaving only a global Z2 symmetry in the resulting theory). Although the 



3 The difficulty in performing the numerical analysis of the TBA equations, required for the 
comparison with CPT, increases rapidly with m. 
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equations have been basically obtained by guesswork and it is desirable to come up with 
further evidence that they are correct, previous experience suggests to take them seriously. 
In the particular case N = 3, where the perturbed Z3-parafermion model is just the ^1,3- 
perturbed 3-state Potts CFT D^~\ the equations for Eo(R) are the same as the ones 
proposed in |21| for the ground state energy in A^~\ corresponding to the flow A$ — ► A4. 

More generally, the ground state energy in Xm does not depend on the MIPF of the 
unperturbed theory; this is predicted by UV-CPT (see sect. 3.2). Assuming the existence 
of the flows within the A-series, it is therefore clear that for all m > 5 Dm^ describes a 
flow of D m to X m -i, but with what MIPF? Our answer is eq. ( |1 . 2| ) . 

In sect. 2 we will describe in more detail known TBA and CPT results, leading to our 
conjecture (|1.2|). In sect. 3 new evidence from UV and IR CPT as well as RG-improved 



CPT is presented. We also discuss in some detail the special case of D$ , arguing that from 
an IR-CPT point of view, i.e. considering it as a ^3^-perturbation (+ higher corrections) 
of A4, it is just the sign of the 4>s,i -coupling that distinguishes it from At- . In sect. 4 we 
briefly discuss our results and methods, concluding with an outlook on some open questions. 
Some observations regarding the appearance of logarithmic terms in CPT expansions are 
described in the appendix. 



2. TBA and CPT results 

We now describe in more detail the results of the TBA approach for the theories 
under consideration. Define the scaling function e(r) = (2tt)~ 1 RE(X, R) corresponding 
to a generic excitation energy E(X, R) in a one-parameter (A, of a specific sign) family 
of perturbed CFTs. Here R is the circumference of the cylinder on which the theory 
is defined (periodic boundary conditions are assumed), and r = MR is a dimensionless 
scaling parameter, M being the mass scale of the theory. M is usually chosen in a way 
that is natural from the point of view of the off-critical theory, e.g. the infinite- volume 
mass of the lightest particle if the the theory is purely massive. It turns out that such a 
choice is also convenient for the TBA equations below. M is related to A through 

|A| = k M y , (2.1) 

where k > is a numerical constant (y is the RG eigenvalue of the perturbing field) . k has 
been determined in many perturbed CFTs numerically, and in a few cases also analytically 
(cf. the appendix). 
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At present, all the TBA-like equations known or conjectured to give the exact scaling 
functions of certain excitations in certain integrable theories are of the "universal" form 

r N roo 

< r ) = "TtE / Mv M H 1 + t ae- £a{6) ) , (2-2) 

a=l J -°° 



where the r-dependent functions e a (9) satisfy the equations 

dff 

6=1 



N f°° c\9' 

e a {0) = ru a (9) -J2 IT K ab (e - 9') ln(l + t b e~^) 



(2-3) 



Here the t a , a = 1, . . . , N (that will be collectively referred to as the "type" t of the TBA 
system) are certain roots of unity, the kernel K is a symmetic matrix whose elements are 
even functions of 9, exponentially decaying as \9\ — > oo, and the v a {9) are of the form 

v a (9) e {m a cosh 9, ±m a e ±e } (2.4) 

where the m a are some non-negative dimensionless parameters (if non-zero, they are mass 
ratios in the theory). 

Ground state scaling functions eo(r) always correspond to the "trivial" type choice 
t = (1,...,1). That other choices of type may yield excitations e(r) in certain theories 
was first proposed in [ 25| |2^] [ 2(| . Fendley |27| then showed that when the theory has some 



global discrete symmetry and is described by a diagonal S'-matrix, the e(r) correspond- 
ing to nontrivial types are the ground state energies of some sectors of the theory with 
"twisted" boundary conditions. These e(r) coincide with scaling functions of excited states 
in the theory (or some "orbifolded" version of it) with periodic boundary conditions. 
With the notable exception of the massless theory A± , the procedure works only for 
the spontaneously broken symmetry phase of the theory, and the excited states obtained 
become degenerate with the ground state in infinite volume. It is not clear at present how 



the approach of [p7fl can be extended to justify the applicability of the "change of type" 
prescription in theories with non-diagonal S'-matrices. 

We will not elaborate here on all the specific choices of N, t a , v a (9), and K a i,{9) 
in (|2.2| )- (|2.3j ) that are known, or conjectured, to give rise to certain scaling functions 
in perturbed CFTs. Rather, we will specialize to the cases we are interested in here, 
namely the theories ■ Let e(Xm \ (p, q)\r) denote the scaling function of the energy 
eigenstate in whose UV limit is created by the spinless primary field <j) p ^ q in the 



CFT X m . This means |2|][§, in particular, that e(xi ±} , (p, q) |0) = djj$ - where 
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d p m q = ^ m 2 "m(m+i) * s ^ ne scann g dimension of (j) p>q and c m is the central charge of 
X m . Recall || that in the CFT A m all the primary fields are spinless and each with 
l<q<p<m— 1 appears exactly once.S The D m and E m models do not contain all of 
the above spinless primary fields. Instead, some primary fields 4>p,q\p,q (with different left 
and right Kac label pairs) with nonzero spin appear, and some spinless fields are doubled 
(the precise field content of the D m models will be given in sect. 3.1). The ground state 
always corresponds to (p,q) = (1,1) = (to — 1, to). 

The equations for the ground state scaling function e(Am \ (1, 1) \r) proposed in [^0 



are 



given by Q-(|J) with N = m - 2, all t a = 1, z^(6>) = cosh (9, u a>1 (9) = 0, and 
K(6) = /( m_2 )/ cosh# where l( N ) is the incidence matrix of the simple Lie algebra A N , 
i.e. = d a ,b-i +8 a ,b+i for a, 6 = 1, . . . , N. The equations for e(Am , (1, l)\r) differ from 
the above only by the choice of i'i(O) and ^ m _2(6>), which are taken to be [^TJ \e Q and \e~ b 
respectively. In [^2| we then conjectured that for m even the equations for e(Am , (2, 2)|r) 
and e(Am \ (^, ^)|r) are obtained from those for e(Am \ (1, l)l r ) simply by changing the 
type to t = (1, . . . , 1, = —1, ts±. = —1, 1, . . . , 1) and t = (—1, . . . , —1), respectively 
(the particular case e(A^ \ (2, 2)\r) of this conjecture appeared also in f2f|). 

In all cases the assignment of the Kac labels to the solutions of the corresponding 
TBA-like equations was justified by an analytic calculation of their r — > limits. Further 
strong support for the conjecture that the solutions give the correct scaling functions was 
provided by numerical calculations and comparison with CPT in the cases m — 4, 5, 6, 7 |22] 
(such analysis of the ground state for m = 4 was first performed in [^Q| ||21|l ). Specifically, 
CPT based on an action of the form ( |1 . 1|) predicts the small r expansion of the scaling 
functions 

oo 

e(r) = e(0) + a n rny + (term(s) nonanalytic in r y , possibly) , (2.5) 

where y is the RG eigenvalue of the perturbing field and the CPT coefficients a n can be 
written as integrated critical correlators (see eq. (3.6) below). In all cases of integrable 
relevant perturbations of nontrivial CFTs studied so far, the possible non-analyticity in r y 
turns out to be given by a single term. This term is either proportional to r 2 , if 2 is not 
an integral multiple of y, or to r 2 lnr otherwise. An r 2 lnr term arises when integrated 



4 Alternatively, we can allow all pairs (p, q) with 1 < p < m — 1 and 1 < q < m, identifying 
(p, q) = (m - p,m + 1 - q). 
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correlators leading to some a n (namely with n = 2/y) still diverge after regularization of 
their UV divergences — required whenever y < 1 — by means of analytic continuation in 
y. Both of these terms can be understood as arising from the RG mixing of the perturbing 
field with the identity operator. The precise coefficients of these terms are determined by a 
renormalization condition, say lim^oo Eq(R) = 0, which is automatically enforced in the 
TBA calculation. No matter what the form of the non-analytic term is, it is clear that it 
is the same for all scaling functions in a given model. It has been determined analytically 
from the TBA equations in many cases. (On the other hand, the part of the small r 
expansion ( |2.5| ) which is regular in r y cannot be derived analytically within the TBA 
approach by present methods. See however, for arguments making its appearance 
plausible starting from the TBA equations.) 

An important consequence of the CPT analysis is that the scaling functions e^\r) 
in the perturbed theory with A = ±|A|, respectively, should be related by analytic con- 
tinuation r y — > —r y , up to the above non-analytic term. This condition, implying the 
relation 

4 +) = (-l) n 4- } (2-6) 

between the corresponding CPT coefficients, already imposes strong constraints on scal- 
ing functions as conjectured within the TBA approach. The ultimate tests are made by 
comparing the explicit values of the (regularized) CPT coefficients, eq. (3.6), with the ones 
obtained from the TBA results. Unfortunately, it is technically very difficult in general 
to carry out the CPT calculations using (3.6) for more than the first one or two leading 
coefficients, due to the lack of convenient representations for the critical correlators in- 
volved;! within the TBA approach, on the other hand, until now it has been possible to 
obtain the expansion coefficients only numerically (to very high accuracy, for the leading 
ones). Though limited, these computations proved sufficient to provide highly nontrivial 
consistency checks on the CPT and TBA results, and as a byproduct accurate values for k, 
of ( |2.1|) were obtained (Table 1 summarizes the results [f20| |21]] ]2^] for the theories relevant 
to this paper; interestingly, the choice of mass scales in Am^ that we implicitly made by 
taking the nonzero y a {0) in the TBA systems to be exactly cosh# and |e ±e , respectively, 
leads to the same K m ||21|| ||22|| in the two theories within the numerical accuracy). 



5 See, however, the generically more powerful method of calculating the a n (numerically) 
based on Hamiltonian CPT, which is discussed in [G2|. 
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III 




Q 
O 


1 

2tt 


4 




5 


0.130234474(2) 


6 


0.11334655(2) 


7 


0.099267(2) 



Table 1: k, eq. (|2.1| ), in the first few models A m . The number in parenthesis is the 
estimated error in the last digit given. (See (A. 9) for a conjectured exact expression for 

K 5 .) 



The analysis of the IR limit r -> oo showed [gIl||22T] that e(A^ h) , (p,p)|oo) = 



(withp 
of the fields d 



= 1 for all m > 4, as well as p = 2, y for even m > 4), exhibiting the flows 
PiP in question in A m to the fields with the same Kac labels in A m _\. These 

More 



flows of fields are expected from the LG and perturbative RG analyses of A 
generally, these methods predict Pfl M for any m > 4 



~rp,P 
-*m — 1 , m — 1 



'p,p 

^2,1 



for p = 1 , . . . , m — 2 



+i for p = 1, . . . , m — 4 



(2-7) 



and in addition the perturbing field ^1,3 is expected JjJ to flow to the irrelevant field ^3,1 
(TT in the m = 4 case [[Tlj) in A m -\. The latter observation lets one hope that it is 
possible to describe the large r behaviour of A^ using IR-CPT, i.e. perturbation the- 
ory around the IR CFT A m _i with the irrelevant perturbation d>^^ (presumably together 
with an infinite series of more irrelevant fields). The large r expansions of scaling func- 
tions obtained using such IR-CPT are expected to be only asymptotic, in contrast to the 
UV-CPT expansions around r = 0, which most probably have a nonzero radius of conver- 
gence [ 17fll|19|l . Nevertheless, successful comparison [^TJ|2^] of the leading terms in these 
asymptotic expansions with the large r behaviour of the conjectured exact results from 
the TBA approach provides further evidence to the correctness of the latter. 

In [^2j we noticed (but at that time did not understand) the following interest- 
ing fact concerning the TBA systems obtained from those of the ground state energy 
e(Am \ (1, l)|f) (see above) with m > 5 odd by just changing the type to t = (—1, . . . , — 1). 
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Denote the corresponding solutions (temporarily) by e^\m\r). Analytic calculation of 
the UV limit gives e^(m\0) = d^+j m+1 — jj. This suggests to identify e^(m\r) as 

2 ' 2 

e(Am j Z2 2^")l r )- However, for the '+'-case this identification turns out to contradict 

the LG prediction of ( |2.7| ), since 

e (+) (m|oo) = cfc'Li-^ 1 • (2-8) 

2 ' 2 1Z 

Furthermore, numerically solving the equations for e^(m\r) for m = 5, 7 we determined 
the small r Vm expansion coefficients a^t\m) (see eq. ( p.5|) ) which we present in Table 2. 
The relation predicted by CPT for a n (A^ } , (^±1, ^)) is not satisfied, though the 
equality (m) = \m) (no minus sign!) is notable. We also noticed in [^2| that in 
absolute value these equal coefficients agree with the CPT prediction (eq. (3.6), where the 
value of K m (see Table 1) that must be used for this comparison was obtained from the 
analysis of the ground state scaling function. 



n 


ai +) (5) 


ai" } (5) 


ai +) (7) 


a ( -\7) 


1 


-0.056643444(3) 


-0.056643449(4) 


-0.042641(1) 


-0.042640(1) 


2 


0.0484328(4) 


-0.0012546(5) 


0.03918(1) 


-0.01743(2) 


3 


-0.020718(4) 


0.01173(1) 


-0.0366(2) 


0.0100(5) 


4 


0.00140(4) 


0.0006(1) 


0.017(2) 


0.008(2) 



Table 2: The first few coefficients in the (regular part of the) small r expansion of the 
functions e^(m|r), m = 5, 7, conjectured to describe scaling functions of the excitations 
specified in eq. (2.9). 



We here propose an "explanation" for the above observations. We conjecture that 
for m odd the correct identification of the above functions e^(m\r) obtained in the TBA 
approach (namely the solutions of (|2.2|)-( p73|) with N = m — 2, t = (—1, . . . , —1), K{6) = 
/(™- 2 )/ cosh 9, and v{6) = (cosh0, 0, . . . , 0) for e^(m\r), v{6) = (±e d , 0, . . . , 0, \e~ e ) for 
e (+)( m | r )) correspond to the following scaling functions: 

e<->(m|r) = e(A&, ^±l)\r) 

(2.9) 

e<+>Mr) = e(Di+\(^±l,^)-\r) . 
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[The superscript '— ' on the Kac label ( 22i 2 ^-, on the second line indicates that the 

corresponding UV conformal state is created by the Z 2 -orfc?spinless primary field d m+i m+i 

2 ' 2 

in D m . This is crucial when m = 1 (mod 4), in which case D m contains two copies of 



* m+1 m-\-l . 
2 ' 2 



In the following we will present evidence for the conjecture (|2.9|) . If true, one concludes 



from ( |2.8| ) and the fact that the perturbing field does not break the Z 2 symmetry, that 
^3 3 02,2 in D { 5 +) : D 5 A 4 



>m + l_ m + 1 ► <Pm-l m-1 1H 

2 ' 2 — o — > — 9 — 



m~\- 1 

2 ' 2 



:D m ^D m . l (m = 7,11,15,...) (2.10) 

2 ' 2 

m-i m-i in Dffl : D m — >■ X m _i (m = 9, 13, 17, . . .) , 



which is our concrete evidence for (|1.2|) when m is odd. (Using IR-CPT we will argue in 
sect. 3.4 that X m _i in the last line of ( |2.10| ) is in fact D m _\. See also sect. 3.3.) Based 
on CPT arguments, we also claim that for the Z2-even spinless primary fields (f> Pjq in the 
model D m , one has 

e(D£\{p,q)\r) = e(A^\ (p, q)\r) (0 M even) . (2.11) 

Similarly, for the Z 2 -odd fields 0m q (q = 1, . . . , ^) in D m with m even 

e(D^,(f,g)-|r) = e{A%\ (f , g)|r) (m = 6, 8, 10, . . .) . (2.12) 

Assuming the LG/perturbative-RG predictions ( |2.7| ) (we do not have any TBA results for 
most of the scaling functions involved) we conclude that the fields in D m corresponding to 
the scaling functions in ( f2.11| )-( p,12| ) flow under the c^i^-perturbation to the fields (possibly 
non-primary) with the same Kac labels as in the A-flows, with the Z 2 symmetry resolving 
possible ambiguities due to doubling of fields in the D-models. To complete the description 
of flows of relevant fields in D m that remain relevant in the IR CFT D m _x, it is natural 
to guess that </>™ +3 m+1 — > <j> m -i m -3 for m = 5, 9, 13, . . ., with the '— ' superscript moved 

2 ' 2 2 ' 2 

over to the IR field when m = 7, 11, 15, . . .. These observations, and the fact that they 



provide support to our main claim ( |1.2|) , will become clearer once we discuss the D-models 
and their perturbations in more detail. 
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3. CPT analysis of the theories Dffl 

3.1. The CFTs D m 

The field content of the models D m is encoded in their MIPFs, written in terms of 
Virasoro characters |HJ || . We have to distinguish between four cases: 

(i) m = 4p + 1 (p > 1): 

4p 2/9-1 2p 

z = J2 J2 ix P , g i 2 +$^ix P ,2p+ii 2 

p=l q odd=l p=l 

2p 4p 2p-l 

+ \x P ,2p+i\ 2 Xp,g X4p+l-p,g • 

p=l p=l q odd=l 

(ii) m = 4p + 2 {p > 1): 

2p-l 4p+2 2p+l 

Z = IXp,g| 2 + ^ IX2p+l, g | 2 

p odd=l 5=1 <?=1 

2p+l 2p-l 4p+2 

+ ^2 l^2p+l, 9 | 2 + 

9=1 p odd=l q=l 



(iii) m = Ap - 1 (p > 2) 



(iv) m = Ap (p > 2) 



4p-2 2p-l 

y y ixp,gi 

p=l g odd=l 

2p-l 4p-2 2p-2 

p=l p=l q even=2 



2p-l 4p 



IXp.gl 



p odd=l q=l 

2p 2p-2 4p 

+ ^^IX2p,q| + Xp,q X4p-p,<; 

g=l p even=2 <?=1 



(3-1) 



(3.2) 



(3.3) 



(3.4) 



This (not always economic) way of writing the MIPFs makes clear the following facts. 
Each term |x p , g | 2 in Z corresponds to (the conformal family of) a spinless primary field 
4> Pt q in the model, whereas a term x v q \t-^ always with (p,q) different from (p, q) 

1 1 



(m — p, m + 1 — q), corresponds to a primary field with nonzero spin. In addition, 

the fields corresponding to the terms on the first (second) line in each case are Z2-even(odd) 
[31]. Note that the models A m and D m (m > 5) have the same Z2-even primary fields. 



On the other hand, the Z2-odd (spinless) primary fields in A m are "replaced" by primary 
fields with nonzero spin in D m , which also contains [^] "extra" odd primary fields which 
are spinless (and lead to a doubling of fields iff m = 1 or 2 (mod 4)). In particular, all the 
relevant primary fields in D m are spinless,! and there are m — 2 such even fields (excluding 
the trivial identity field </>i,i) and two odd ones; in A m , on the other hand, there are 
2(m — 2) nontrivial relevant (spinless, primary) fields, half of them even and half odd. 
An important ingredient in the CPT analysis below is the structure of the operator 



algebra of the models D m . The fusion rules for these models were given explicitly in pi 
They are essentially the A-series fusion rules (applied separately to the left and right Kac 
label pairs (p, q) and (p,q) appearing in (|3J])-(p\l])) intersected with the Z2-symmetry 
selection rules. On the other hand, the problem of finding the (nonvanishing) operator 
product expansion coefficients (OPECs) is much more involved ]32|] than in the corre- 



sponding A-models. This is the case for CFTs with non-diagonal MIPFs in general [33 
The main subtlety lies in the signs of certain OPECs in the non-diagonal models. What is 
clear [33 [[33] , however, is that if the fusion of two spinless primary fields in the non-diagonal 



model gives rise to primary fields that are all spinless, then the OPECs involved are iden- 
tical to those in the diagonal model.@ For D m , it is known that the nonvanishing 
OPECs (of primary fields of definite Z2 parity) satisfy 

(Pl,9i;Pl,9l)(P2,92;P2i92)(P3,93;P3!93)/ (Pl,<?l)(P2,<22)(P3,<?3) (pi , 9i ) (p 2 , Q2 ) (P3 > 93 ) ' ^ * ' 



6 A notable exception is the 3-state Potts model D$ that contains two relevant (d = |) primary 
fields with nonzero spin: 02,i ; 3,i and <^3,i ; 2,i- 



More generally, if the correlator of some fields in conformal families whose ancestors are all 
spinless primary fields gets contributions from intermediate fields (through "factorization") all 
of whose ancestors are again spinless, then this correlator is the same as in the corresponding 
CFT with diagonal MIPF. This follows from the uniqueness of the solution to the monodromy 
problem p3| ]. 
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where the OPECs on the right-hand side are those of the primary fields in A m |3|[ . 



Now recall that the primary fields in A m are normalized, up to signs, by requiring that 

((t> Pl>qi (00,00) <^>p 2 ,g 2 (0,0)) = 5 PlP2 5 qiq2 (here 1 < q { < p { < m - 1). The signs can be 



further fixed so that all the OPECs in A m are non-negative [ 34j] . Hence eq. (|3.5|) gives us 



the OPECs in D m up to signs, which turn out to be crucial for the discussion below. 
3.2. UV-CPT for 

We are now ready to use CPT to study the theories Dm\ trying to support the 
validity of eqs. ( |2.9| )-( p7T2|) that are our evidence for the main claim ( |1 . 2|) . For the scaling 
functions e{X^\ (p, q)\r), CPT based on ([□) gives the small r expansion (|2.5|) with the 
CPT coefficients 

a n (x£\(p,q)) = -(2tt) 1 -^ ( TKn ) n 



n — 1 ,2 n — 1 



X 



J II (2%\z-\)ym (^g( 00 ' 00 ) 01,3(1.1) J[ ^(zj.Zj) Pi ,(O,O)) c 



(3.6) 



The correlators here are connected (with respect to the "in- and out-states" created by 
4>p,q) ( n + 2) -point functions in the CFT X m on the plane. [On the r.h.s. we suppressed 
the factor ((j) Pjq (oc, 00) <p Pjq (0, 0)) _1 , which is set to 1. In particular, we will assume this 
standard CFT normalization for the combinations ^ q of definite Z2-parity in the relevant 
.D-models. Since in the perturbations of these latter models (like in all the other 
and D^' theories) the Z2 charge is a good quantum number also away from criticality, 
it makes sense to consider states in the perturbed theory whose UV limits correspond to 
(f) p ^ q and 4> Pjq ] The problem of (UV and IR) regularization of such integrated correlators 
was discussed at length in |22|] (see also the appendix). 

Of particular importance to us will be the first CPT coefficient, which is never diver- 
gent, given simply in terms of conformal OPECs by 

aiP^,^)) = ±(2tt) 1 -^ Km G%$ lt3)M ■ (3.7) 

Note that it always vanishes for the ground state (in the unitary theories under con- 
sideration) corresponding to (p, q) = (1,1). 



Let us now prove eq. ( ^.11| ). This equation is a consequence of the fact that all spinfull 



primary fields in the D-models are Z2-odd and therefore cannot appear in the fusion of 



13 



two Z2-even fields. It is therefore clear (cf. footnote 7) that the CPT coefficients involved 
are identical, including signs, in Dm and Am . Analytic continuation of the small r CPT 
expansion to all r proves ( |2.11|) . Actually we see from the above argument, which did not 
use the fact that the Z 2 -even fields <p Piq involved are primary, that eq. ( |2.11|) generalizes 
to the statement that the whole Z2-even sectors of the finite-volume spectra of D^ and 
are identical. 

Similarly, ( |2.12| ) follows from the fact that in the fusion (in D m , m even) 

[01,3] x = [</>^ 9 _ 2 ] + + [<% g+2 ] (3.8) 

no primary fields of nonzero spin can appear on the r.h.s. (there are no such fields for which 
the first entries of the left and right Kac label pairs are both m -). Again, this implies that 
the relevant CPT coefficients in Am^ and Dm are identical (this is true also for all the 
scaling functions corresponding to descendants of the primary fields ^ q ); the fact that 
(f>m q is Z2-even in A m if m = 2 (mod 4) does not matter. Note that in this latter case 
( |2.11| )-( p.i2D imply, in particular, that the Z 2 -doublet ™ in D m flows to the doublet 
with the same Kac label in D m -\. 

Now to ( |2 . 9| ) — ( pTT0|) , our most interesting claim. Note first of all that in the conjectured 
flow (b m+i m+i — > <j) m-i m-i (with appropriate ' — ' superscripts as in (E.IOI)) in Dm — > 

2 ' 2 2 ' 2 1 1 

D m _i, m odd, the field becomes more relevant during the flow, i.e. its scaling dimension 
decreases: 

d {m - 1] - d {m) - -— m2 + 3 < (3 9) 

" m — 1 m-1 u m + 1 m + 1 . n 1 \' J '' J J 

2 ' 2 2 » 2 4771 771 — 1 

This is in contrast to all the known flows of fields in the A-series, where fields become less 
relevant (the fact that this feature of the A-series flows, reminiscent of the decrease of the 
central charge along unitary RG trajectories |TJ, is not universal, was already noted in JTT|). 
The qualitative difference can be seen already in the first order of the small r expansion 
of the relevant scaling function: Remember our empirical observation concerning the first 
expansion coefficients of the "scaling functions" e^(m\r) of the type t = (— 1, . . . , — 1) 
TBA systems with m odd. Assuming ( |2.9| ) we can now state it equivalently as 



+ ai(A(-\(^^)) 



a i 



(Am \ ( 



m-\-\ m-\-\ 



(m odd) 



(3.10) 
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where the '— ' superscript on the l.h.s. can be ignored in the case m = 3 (mod 4). Actually 
we found [^2] the first equality (the second is just a particular case of Q2.6p) only for 
m = 5, 7, but in the meantime also have checked it for m = 9, 11, and believe that it holds 
for all odd m. 

What is the origin of this difference between the flows in the D and the A series from 
the point of view of CPT? We claim that for m = 3 (mod 4) the crucial fact, proving 

(fng), is 

Q-V m + 1 m+lu, ,wm + l m + 1 \ = ^/ m+l m + lvi n\ / m + 1 m+li 7 171 = 3 (mod 4) . (3.11) 

\— — ' 2 — J^WA — 2 — ' — 2 — ) ( — 2 — ' — 2 — A 1 ) 3 ^ — 2 — ' — 2 — > 

The crucial sign difference is consistent with the fact that now spinfull fields appear in 
the fusion [01,3] x [0 m+i m+i ].l Note that sign differences between OPECs like the one 
in ( |3.11|) have a less trivial effect on higher CPT coefficients, so that for n > 1 we do not 
expect any simple relations like ( |3.10| ). This is consistent with our numerical results — 
cf. Table 2. 

For m = 1 (mod 4), on the other hand, the crucial fact is, we claim, 

(T < -~ D " 1 * ) — _(n( D m) — 1 ( rr . nr \ A\ 

^ ( m + 1 m + 1 \ _ n qVHi+l m + 1 \— ^ I m+1 m + 1 \+fi o\ ( m + 1 m + 1 ^ 5 " b ± L ± ) 

\ 2 ' 2 > \*"> a l\ 2 ' 2 I \ 2 ' 2 I \*"> a I\ 2 > 2 ' 

_ _(p{A m ) 

V 2 ' 2 2 ' 2 / 

This allows the scaling dimension to decrease in the flow ~ +1 m+1 — > m — l T7t — l in Drn ^ , 

— ' —2" 2 ' 2 

whereas in the flow 0~t .-, „ — >• m-j-1 m-f 1 the scaling dimension increases, as in the 

— 2~ >- 2~~ 2 ' 2 

corresponding A-series flow. 

For m = 5 we can use the Z3 symmetry of the 3-state Potts model -D5 to easily prove 
the first equality in ( |3.12| ) (note that the second equality holds in general, according to 
our remarks in sect. 3.1). Under this Z3-symmetry 03 ; 3 = -^(0^3 + ^3,3) an d 03,3 = 
-^(03^3 — i<Ps 3) form a doublet of oppositely charged fields. Recall that the perturbation 
is by 0j 3 = ^71 (0i,3 + 0i, 3 )) which implies 

/p(-De) _ _L^_rp( D 5) _ _A_J_tf( D s) fqiq^ 

^(3,3)±(l,3) + (3,3)± - =C 2^^03,3±^,3' 01,3+01,3' 03,3±0J,3 ~~ yM 01,3, 03,3 ' ^ ^ 



Unless we misunderstand the results of [p2| , our ( 3.11 ) disagrees with them when 



m 



3 (mod 8). The same remark applies to our (3.12), now with m = 1 (mod 8). 
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where we used the linearity of the OPECs and in the last step also the conservation of 
the Z3 charge and the Z2 symmetry (we choose to assign the same Z3 charge to and 
^1,3). This proves (|3.12|) for m = 5. 



[It is interesting to explicitly see how different the consequences of the Z3 symmetry 
are in the theory Dg . (Dq, the tricritical 3-state Potts CFT ||, and D$ are the only Z3- 
symmetric models among the X m .) First, note that there is a single Z3- (and Z2-)neutral 
field (pi y 3 in Dq, in contrast to -D5, so that both the Z2 and the Z3 symmetries are preserved 
in the perturbed theories (This fact by itself shows that cannot possibly flow 

to A$, as A§ is not Z3-symmetric.) Now consider the first CPT coefficients relevant for the 
flow of the doublet (p^ 3 — > 4>f 3 in Dq + \ Here, in contradistinction to ( |3.13| ), conservation 
of the Z3 charge gives 

(ni D 6) _ (n( D e) _ (n(A>) (a 1 a\ 

^(3,3)±(1,3)(3,3)± ~ t J%±*; iS , 4>i,3, *s,s±0; >a ~~ %3,3, tfi.a, ' K^'l-V 

implying that ai(.Dg + \ (3, 3) + ) = ai(.Dg + \ (3, 3)~) as expected from other considerations 
discussed earlier.] 

3.3. RG-improved CPT 

To gain direct information about the IR fixed point from standard UV-CPT on the 
cylinder, as discussed in the last subsection, would require analytic continuation of the 
small r expansion to large r. This is generically not possible in practice. However, one can 
sometimes directly "see" the IR fixed point by using the RG to (partially) sum up pertur- 
bation theory. This assumes, of course, that one knows the /3-function quite accurately up 
to its IR zero. In the present context this restricts the applicability of this approach to 
large m, where the IR fixed point is close to the UV one (as measured by the renormalized 
coupling, not the bare coupling A of sect. 3.2 in terms of which the IR fixed point is at 
A = 00). This method has been discussed in detail in the literature P| P] [J13| , so there is 
no need to review it. We just summarize the results relevant for us. 

Consider the RG flow induced by perturbing a CFT of central charge c by a slightly 
relevant (y 1) operator <fi, which up to irrelevant operators closes on itself and the 
identity operator under repeated fusions. The central charge d at the new IR fixed point 
(attained by choosing the appropriate sign for the perturbing term) then satisfies 

c'-c = -^A? + °(y 4 ) ■ (s- 15 ) 
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For a field </> a which to leading order does not mix with any other field (like the fields PiP 
on the "diagonal" of the Kac table of a minimal model) the change in scaling dimension is 

d' a -d a = ^t^y + 0(y 2 ) . (3.16) 



In the case we are interested in </> = </>i j3 , y = = -^= + 0(±) so that 

d — c = — ¥5 + O(^z), consistent with a flow X m — > X m -\. For </> Q = cf) p>p with p < [^] 

2 — 1 2 im 

odd the OPECs G^ a ^ a = 2 P ^ m2 + O(^g) are identical in A m and D m , showing the 
existence of flows <^> P)P — > (/> PiP in X m — > X m -\. These flows do not allow us to directly 
distinguish between, say, D m — > D m -i and D m — > A m -i (although it is clearly hard to 
see where the fields </> PjP , p even, in A m -i should come from!). 

For m = 3 (mod 4), however, the flow of the Z2-odd field 4> a does enable 

2 ' 2 

us to directly distinguish the two possibilities for flows starting from D m . Depending on 
the sign of = ±(^75 + ®(^)) we now have 

d! a -d a = ±± + O(^), (3.17) 
indicating the flow of scaling dimensions d { ™+i m+i —> d ^+^m+i or d ( Z+i m±i 

2 ' 2~ ~ 2 ' 2 2 ' 2 

— > d ^Zj_ m=i , respectively. According to our claim (|3.11| ) the second possibility holds 

for Dm , which is only possible if D m — > -D m _i (because of the Z 2 symmetry, as there is 

a single field (bm-\ m -i — which is Z2-even — in A m -i). 

2 ' 2 

For m = 1 (mod 4) the analogous sign difference (|3.12|) again indicates the flow 

(P m + l m+1 > (J) m-1 m-1 lli D£\ But now, unfortunately, this flow alone does not allow us 

2 ' 2 2 ' 2 

to decisively conclude that the IR CFT is _D m _i. In fact, for m = 5 we know this cannot 
be the case. To show that D m — > A m -\ is not possible for m = 9, 13, 17, . . . requires the 
study of flows of other fields, not on the "diagonal" of the Kac table, but we will not 
pursue such studies here. Intuitively, the fact that there are about half as many relevant 
fields in D m than in A m _i, for m ^> 1, makes the flow between these two theories highly 
implausible. A more explicit version of this argument will be offered in sect. 3.4. 

For m even the OPEC (C^ a( ^ Q , 4> a = ^^,^5 i n D m (with superscripts '±' for m = 
2 (mod 4)) is identical to the corresponding one in A m . Therefore (f)m m in D m , m 

2 (mod 4), must flow to fields 6% m in X m _i, which is only possible if X m _i = D m _i. 

2 ' 2 

For m = (mod 4) there is only one (Z2-odd) field 0>n ) >™., and so this argument cannot be 
applied. As in the case of m = 1 (mod 4), it is necessary to study the flow of other fields 
to conclude within the RG-improved CPT approach that D m — > is not possible. 



17 



34- IR-CPT for 

We finally discuss what can be learned about the RG flows we are considering from 
CPT around the IR fixed point. The basic idea is that the action 

A Am _ x + J d 2 x [g 03,i (a;) + ... ] , m>5 (3.18) 

can be used to study perturbatively the IR asymptotics of the theory Am\ in particu- 



lar |21 1 [2^] the large r behaviour of the energy scaling functions in this theory. The leading 



nontrivial term in these scaling functions will generically be proportional to r~ , since 
d^l = 2 + The '. . .' in the integrand refers to a presumably infinite series of 

fields more irrelevant than 03,1. (Note that the perturbation theory we are considering 
here corresponds to a non-renormalizable interaction in standard Lagrangian QFT.) Given 
the first (or perhaps more, see below) terms in ( |3.18| ), the infinite and finite parts of the 
counterterms written as '. . .' are determined, in principle, by two conditons: a) That they 
make Eq(R), say, finite order by order in perturbation theory, and b) that the presumably 
asymptotic large r series is consistent with the TBA and/or the (analytic continuation of 



the small r) UV-CPT series. In [^TJ |2!| it was shown that the first few terms (independent 
of '. . .') in the large r expansion of certain scaled energies for m = 4,5,6,7, calculated 
from (|3.18| ) as in UV-CPT, agree with those obtained in the TBA approach. 

We will now examine what ( |3.18| ) and its generalization to the .D-flows can teach 
us. For m > 6 we should consider, according to ( |1 . 2| ) and ( |2.11| ), the action ( |3.18[ ) with 



A m -\ replaced by -D m _i (and <p^,\ taken to be 4>^ x when m = 7; note that in this case 
the IR perturbation breaks the Z3 symmetry of Dq, which explains the accidental - 
from the point of view of UV-CPT — appearance of this symmetry at the IR limit of the 
theory D^ +) )E The comparison between the "backward flows" induced in the D-series and 
those in the A-series essentially parallels the previous discussion of UV-CPT. Again, the 
conclusions are that flows of certain fields are the same in the two series (provided the IR 
perturbation is identical), whereas for other flows (e.g. (|2.10|) ) the difference is consistent 
with differences in the operator algebras of the models -D m _i and A m -\. 



9 A similar situation occurs in the flows of A5 and D5 to A4, the (superconformal) tricritical 
Ising CFT. There the IR perturbation (f>3,i, the bottom rather than the top component of (the 
determinant of) the super-stress-energy tensor, explicitly breaks the supersymmetry. 
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The most interesting case is m = 5, on which we now elaborate. Here IR-CPT seems 
to present us with a little puzzle: If we are correct in our claim that is the common 
IR limit of both A { 5 +) and D { 5 + \ and moreover that the IR fixed point is approached in 
both cases along the direction of 03, 1, what is the difference in the IR perturbations of A4 
that leads to the two different "backward flows" to A 5 and D5? Recall that by UV-CPT 
we concluded that the Z2-even sector of the spectrum of the two theories is identical, so 
that the IR limits of UV fields in the conformal families of p ,i (p = 1,2,3,4) and 0^~ 3 
(p = 1, 3) in D§ and A§ (the superscript '+' being redundant in the latter) is the same. In 
particular, by (|2.7f) we have for the relevant Z2-even UV fields 0i,i — > 0i,i, (p 33 — > 03,3, 
02, 1 - > 01,3) an d 0^3 — > 03,i- However, the Z 2 -odd sector is necessarily different: Our 
main TBA result ( |2.9| ) implies that 0J 3 — > 02,2 in ^5 whereas ( |2.7|) predicts 02,2 — * 02,2 
in , ( |2.7| ) also predicts 0x,2 — > 02, 1 in A^~\ and our guess for the "corresponding" flow 
in is 0^~3 — > 02,i- We also believe that all the other Z 2 -odd fields in D 5 , including 

the primary spinfull ones, flow to descendants in the families [02,2] and [02, 1] in A4, the 
UV origin of these families in A 5 being [<j> P} 2[ (p = 1, 2, 3, 4). 

Trying to understand the difference in the IR perturbations leading to D$ and , 
the first thought that comes to mind is to blame the unspecified '. . .' for this difference. 
However, this would imply that the 0{r~ l ) corrections to the IR limits of all the energy 
scaling functions in the two theories are the same, whereas based on the flows described 
above one would expect this only for the Z 2 -even sector! We would like to suggest that 
actually the IR perturbations are already different in the first 03,i term, the difference 
being the sign of the coupling g. [Note that this sign difference is also the only way to 
avoid the following potential contradiction: If the leading IR perturbation were the same, 
then conditions a) and b) mentioned earlier and the equality of Eq(R) in A^ and 
would imply that these are described by exactly the same IR perturbation of A4, i.e. are 
the same theories, which is wrong. There is however one caveat, namely, if ( |3. 18|) contains 



is 



terms of the form g a J d 2 x(f)(x) : a > 1, which are distinct for A 5 and , and 
such that it does not appear in the (repeated) fusion of 03 ; i with itself. We consider the 
existence of such terms unlikely, but do not really have an argument to exclude them.] 

Consider, in fact, the perturbative expansion of scaling functions based on ( ^.18| ) for 
m = 5 with the '. . .' ignored, namely Aa a + g J 03, i- The resulting CPT coefficients are 
of the form ( |3.6| ) with 0i ; 3 replaced by 03, 1, y m by 2 — d 3 \ = — 1, and K m by K m . Now 
as already noted in |2~2" |, in A4 03 i is the same as 0i,4, the latter being a Z2-odd member 
of the (1, (^-operator subalgebra of the model. (This Z2 symmetry, corresponding to the 
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self-duality of the tricritical Ising lattice model, has nothing to do with the Z2 symmetry 
discussed so far in the paper, which corresponds to "spin reversal". In particular, the 
(1, (j)-operator subalgebra of At, containing both Z 2 -even and Z 2 -odd fields, constitutes 
the whole Z2-even sector of the model.) Therefore all correlators (<fi Hi 01,4 0) with 
a Z2-even field 4> (> e - °f definite Z2-charge) and n odd vanish in A4, and so do the 
corresponding CPT coefficients. As a result, the whole perturbative expansion of scaling 
functions corresponding to the Z2-even IR fields is independent of the sign of g, and is 
in fact in powers of r~ 2 (modulo an r _2 lnr term, see below). On the other hand, for 
the Z2-odd IR fields in the conformal families of 4*2,2 and <p2,i (not belonging to the (1, q) 
subalgebra, i. e. not being Z2 eigenstates) the sign of g does matter. 

Of course the above arguments have to be generalized to include additional IR perturb- 
ing fields. [One may try to exploit our observations regarding the finite- volume spectra of 
the two theories X§ when trying to determine these additional perturbations. For exam- 
ple, allowed perturbing fields are all the spinless descendants of </>3,i, whose RG eigenvalues 
are necessarily odd negative integers, and all spinless descendants of the Z 2 -even identity 
field, whose RG eigenvalues are even negative integers.] Still, they are very suggestive and 
make us believe that in the (conjectured) RG flows from A5 and -D5 to the common IR 
fixed point A4, the latter is approached from exactly opposite directions. 

An important lesson of the above discussion is the following. For m > 6 the n-point 
functions of the IR-perturbing field in X m _i do not vanish for all odd n. Hence, 
consistency with the UV-CPT observation that the Z2-even sectors of Am and Dj +; 



' rn 



have the same finite- volume spectrum, seems to require that the IR perturbation of AT m _i 
leading to Am^ and Dm is the same, also in sign. Therefore, Am^ '■ A m — > A m -i and 
Am 7^ clearly imply that D m flows to D m -i, not A m -±, for m > 6. By the same 

argument, the CFT at the IR limit of Em \ m = 12, 18, 30, must be E m -\. 

To conclude this section we briefly describe results of the TBA approach relevant to 
the IR behaviour of . In this case, the energy scaling functions for which integral equa- 
tions have been conjectured are e(A^~\ (1, l)|r) = e(D^ + \ (1, 1) |r) and e(D^ + \ (3,3)~|r). 
Solving numerically the TBA equations, we found [^2| for r ^> 1 

e(A^ +) ,(l,l)|r) = -— - 0.02723(2) r" 2 lnr + 0.0173(2) r" 2 + . . . (3.19) 

for the ground state, and for the excitation describing the flow <p^ 3 — > 02, 2 m we 
estimate 

e(l4 +) ,(3,3)-|r) = -L + 0.03333(1) r" 1 + .. . . (3.20) 
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These fits are based on numerical results for the scaling functions evaluated at 200 < r < 
600. 

We should emphasize that in general it is much more difficult to perform the fits of 
the TBA scaling functions at large r, leading in our case to (|3.19|) -( |3.20| ), than those of 
the small r dependence leading to UV-CPT coefficients (Table 2, for instance). We were 
therefore able to determine with some confidence only the leading (asymptotic) expansion 
terms indicated in (|3.19|) -( |3.20| ). Still, the basic features predicted by IR-CPT based on 
( |3.18|) are already noticable, namely the expansion in powers of r -1 . The presence of the 
r -2 lnr term in ( |3.19| ) indicates the divergence of the second CPT coefficient (multiplying 
r -2 ) even after its UV regularization via analytic continuation in y m = 2 — ^31 ■ As in UV- 
CPT, the same non-analytic in r mn term is expected to appear in the large r expansions of 
all the scaling functions in the model, in particular in ( |3.20|) . Unfortunately, the accuracy 
of our numerical results for e{D^~\ (3, 3)~|r ^> 1) does not allow us to test this prediction. 

There is however one intriguing observation we can make based on (|3.19| )-( j3.20|) , which 
can be taken (with a grain of salt) as a consistency check of these results. First, using ( |3.7| ) 
with C(^2)(3 i)(2 2) = M 10 we ex ^ r act from the C(r _1 ) term in ( |3.20|) the value k ir = 
14n~ 2 ■ 0.03333(1) = 0.04728(2). We now extend the observations made in the appendix 
regarding logarithmic terms in UV-CPT to the case at hand, using the recipe (A. 4) to 
treat the logarithmic divergence of the leading IR-CPT expansion term of the ground 
state. The divergence is now due to the simple pole at y = — 1 in eq. (A.l), whose residue 
is Rq = — 7t 4 k 2 r /16. Eq. (A. 4) then expresses the coefficient of the r -2 lnr term in ( |3.19j ) 
as B in = aRo, and from the numerically obtained values of K in and B = —0.02723(2) we 
compute a = 2.001(3). After reading the appendix, the reader will hopefully believe that 
in fact a = 2, which leads to the (conjectured) exact relation 

B IR = . (3.21) 



4. Discussion 

We have presented evidence for the existence of RG flows between members of the 
.D-series of minimal unitary CFTs. These flows are induced by perturbation of the UV 
CFT by the least relevant (spinless) field in the model, a perturbation that preserves the 
integrability of the theory away from criticality. We have studied the explicit flow of scaling 
dimensions from the UV to the IR CFT for various operators. For the RG flows D5 — ► A4 
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and A$ — > A4 we have presented strong evidence that they approach A4 from exactly 
opposite directions. 

If there is one general lesson to be learned from this and other related studies, we 
think it is this: In trying to understand integrable QFTs interpolating between different 
RG fixed points, the study of the finite-volume spectrum can be a very powerful tool. The 
small- and large-volume behaviour of the spectrum gives information about the UV and 
IR fixed points. It can be studied with a variety of analytical and numerical techniques, 
like the TBA, UV- and IR-CPT (as well as the "truncated conformal space approach" |36| 
which we have not utilized here). Admittedly, "TBA" integral equations (in particular 
for excited states) often involve some guesswork, but the structural regularities observed 
in the known TBA equations make such guesses quite natural, in many cases. The TBA 
approach is non-perturbative and the equations for excited states show quantitatively how 
certain fields flow from the UV all the way to the IR; we think this is at least as interesting 
as the LG analyses (which in many cases are only qualitative, at best) and the perturbative 
RG calculations which only apply if the theory is "close" to a model that can be expressed 
in terms of (compactified, perhaps) free fields. 

We would like to comment on an interesting consequence of ( |2.9| ) regarding the "stair- 
case model" of Al. Zamolodchikov. This model is given by a one-parameter (#0) family 
of diagonal S'-matrix theories (the status of the underlying QFTs is not clear) describing 
the scattering of a single particle. The TBA analysis of the model shows |35| that as the 
parameter 9q gets larger, the ground state energy scaling function approximates better and 
better all the functions e(A+, (1, l)|r), m = 3, 4, 5, . . ., "stringed together" (see |35| for the 
precise meaning of this statement). Therefore in the limit 9q — > 00 the model describes, in 
some sense, the whole RG trajectory of flows between the unitary minimal models starting 
at m = 00 (c uv = 1) and going all the way down to the Ising CFT m = 3 and finally to 
the trivial massive IR fixed point (c IR = 0). 

In |27j the TBA equations for the lowest energy level of the model in the sector of anti- 



periodic boundary conditions were considered. Using the results of J22[ the corresponding 
scaling function was seen to approximate (as 9q — > 00) the whole set of e(Am , (^rS ^)l r ) 
for m even and e^ + \m\r) for m odd, stringed together (in an alternating pattern). As was 
pointed out in |Z7| , this does not approximate a sequence of excitations in the A-series flow. 
But now having the new insight ( |2.9| ), our interpretation is that the sequence of excitations 
involved is in the .D-series, describing the flow of spin fields . . . — > 4>2p+i,2p+i ~* fap^p ~^ 

4>2p,2p -> 02p-l,2p-l 02p-l,2p-l -»• ■ • • ^3,3 ~+ 02,2 -> 02,2 in • • • ~> D 4p+1 -> 
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D^p — > D^p-x — *■ D/±p_2 D^p-s — > . . . — > D 5 — > A 4 — > A3, respectively, and finally in 
the last step the relevant excitation is e(A 3 , (2, 2) |r) that becomes degenerate with the 
ground state e(A$ , (1, l)\r) = e(A 3 (1, 1) |r) in the IR limit r — > 00. 

The significance of this observation is not yet clear. In particular, does it indicate 
that the staircase model actually interpolates (asymptotically) between all the CFTs in 
the .D-series rather than in the A-series (recall that the finite-volume ground state energy 
is the same in the A- and D-flows)? Or, does it indicate the existence of an "orbifolded 
staircase model" approximating the .D-flow, with the modified TBA equations discussed 
above describing an excitation in this model? Obviously, to answer these questions a better 
understanding of the staircase model itself is required. 

Related to this, the reader might wonder why the simplest possible integral equation 
for an excited state in \ namely the type t = (—1, . . . , —1) system discussed in sect. 2, 
describes an excitation in Dm and not Am^ when m is odd (for m even the excitation 
energies in question coincide in the theories A^ and Dm)- We do not have a satisfactory 
explanation of this at present, it is basically an "empirical" observation. However, we 
note that in the study [28] of off-critical orbifolds of the lattice models related to Am a 
fundamental difference was observed between the cases of m even and m o dd0 (specifically, 
applying the "orbifolding" procedure |28] to the lattice models of Am leads to a model in 
the -D-series, namely -D^+ 3 , only when m > 5 is odd). We think that this work, together 
with P7| , might be helpful in understanding the issue. 



Note added 

While putting the final touches on this paper we received |3!J, which presents a per- 
turbative argument for flows between minimal models high up (m ^> 1) in the .D-series. It 
is based on [ |D[ , where the difference of the torus partition functions of Am at its UV and 
IR fixed points is calculated to leading order in 1/m. The fact that flows in the D-series 
can be distinguished from those in the A-series already to leading order, turns out to follow 
essentially from the fact (cf. our sect. 3) that there are half as many spinless primary fields 
in D m than in A m , for large m. To leading order, at least, it is not possible to follow the 
flow of individual operators using this method. 
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Appendix A. Logarithmic terms in CPT 

We present some observations regarding logarithmic terms in the small r expansion 
of energy scaling functions in certain perturbed CFTs. As mentioned in sect. 2, the 
appearance of such terms is indicated within the framework of CPT by the divergence, even 
after analytic continuation in y (the RG eigenvalue of the perturbing field) , of certain CPT 
coefficients in the perturbative expansion in powers of r y . It seems that CPT is completely 
useless in such cases, in particular there is no way to predict the form of the resulting non- 
analyticity in r y in the true non-perturbative answer. However, by considering cases where 
exact non-perturbative information is available, we empirically arrive at a recipe which - 
though incomplete, as we shall see — allows one to obtain some interesting (conjectured) 
results. 

Consider a generic perturbed (unitary) CFT defined by an action of the form (|1.1|), 
with a single perturbing field <j> of RG eigenvalue < y < 2 replacing <pi^ there. (We 
assume that (p is the most relevant field — except for the identity — in some subalgebra 
of the operator algebra of the model, so that no nontrivial fields have to be added [0 as 
counterterms to the perturbation when renormalizing the theory.) For the ground state 
energy scaling function eo(r) one can obtain the two leading CPT coefficients analytically 
from (|3.6| ) (see e.g. P2| ). Namely, in the CPT expansion eo(r) = — ^ + XmLi a nf ny we 
have a i = (cf. ( |3.7|) ) and 

a 2 (y) = -I(27r) 2 ( 1 -^K 2 7 2 (l-!)7(2/-l) , (A.l) 

aa(») = C w7 3 (V)7(¥) ■ ( A - 2 ) 

Here ■y(s) = r(s)/r(l — s) and C^^^ is an operator product coefficient. (In order not to 
complicate the notation we considered only the perturbation with positive A.) 

Eqs. ( |A.lp and ( JA.2[ ) actually give the analytic continuation in y of the correspond- 



ing integrals (|3.6|), which converge only when 1 < y < 2 and | < y < 2, respectively. 
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Known non-perturbative results of the TBA approach [2C][21][22] lead one to conclude 
that these analytic continuations (and analogous ones for higher CPT coefficients that we 
are unable to compute explicitly) provide a consistent renormalization scheme — not only 
UV regularization — for all y E (0,2) different from jt, where iV is a positive integer. 
(We are not aware of a rigorous proof of this fact, though it is highly plausible, since this 
kind of renormalization scheme is essentially equivalent to dimensional regularization with 
minimal subtraction in ordinary QFT.) 

But what if y = jt, N = 2, 3, . . .? Then the analytically continued a N has a simple 
pole at jt (a n with n < N are finite after analytic continuation and those with n > N 
are convergent to begin with — see Jl9j [ 22| and references therein for detailed discussion) . 



The simple pole indicates that if one introduces a UV cutoff e to regularize the divergence, 
then the integral for a N diverges logarithmically when e — > 0, unlike powerlike divergences 
in e that are encountered in a n with n < N. There are two cases, the theories A^ and 
A^ where V = jt with iV = 2, 3, respectively, in which in order to obtain the known 
exact TBA results the following recipe can be used: Expand (the analytically continued) 
a N in Laurent series around j? 

oo 

a N ( y ) = (y~^ k ; ( A - 3 ) 

fe=-l 

then in the CPT expansion of eo(r) replace 

a N r Ny = a N r 2 with (aR ln/3r + i^) r 2 (A.4) 

where a, (3 are certain (real) constants. 

The remarkable fact is that a comes out to be a "nice" rational number — see below.0 
Moreover, the recipe is supposed to be "universal" in the sense that one has to use the 



11 This is also the case in a similar recipe that was proposed by Dotsenko f?7f for treating 
divergences in the CPT expansion of the spin-spin correlation function in the so-called Ising field 
theory (IFT) on the plane. There the scaling variable analogous to our r is proportional 

to the separation of the fields rather than the volume of space. A major difference between 
Dotsenko's problem and ours is that his recipe is meant to cure IR divergences that are present 
in infinitely many orders of CPT, whereas we deal here with UV divergences of only finitely many 



terms in the perturbative expansion. This is manifested in the non-perturbative result [38] for the 



problem in [37] by the existence of infinitely many power-log terms in the expansion, whereas in 



the finite-volume energies there is (apparently) at most a single such term. 
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same a and f3 when treating the logarithmically divergent iV-th CPT coefhcents for all the 
energy scaling functions in the given model. This requires that we keep the R± term in 
( |A.4|) , which could have been absorbed in (3 as long as only the ground state is concerned, 



allowing for different coefficients of r 2 in different scaling functions (note, though, that the 
coefficient of the r 2 In r term is the same in all of them, as follows from the fact that the 
"strength" of a possible UV divergence in ( |3.6|) is independent of the in- and out-fields 

We first demonstrate the use of our recipe forthelFT A^K In this case the complete 



finite- volume spectrum is known analytically (see sect. 6 of [IIT1 and references therein). 
For the ground state scaling function, in particular, 

e (r) = - — r 2 lnr + — In (ne^~ lE \ r 2 + C(r 4 ) , (A.5) 

24 2 2 V / 

where "f E = 0.577215 ... is Euler's constant and k = (k is exactly known from the 

Lagrangian formulation of Ag as obtained from a free massive Majorana fermion through 

(+) 



a "GSO projection" P^j , since the perturbation Xcfti^ leading to A 3 is equivalent there to 



the fermion mass term Y-itptp, in conventional complex notation.) Using this exact result 



we conclude from ( |A.1|) and (|A.4|) that 



a = 2 , f3 = i e lE ~^ (A.6) 

in IFT. We verified that our recipe gives the correct r 2 term also in the scaling function cor- 
responding to the field TT of the UV Ising CFT (the integration of the relevant correlator 
can be performed analytically using results of sect. 4.2 of ET 



We turn to the ground state in Ag , for which Al. Zamolodchikov evaluated the coef- 



ficient of the r 2 In r term analytically from the TB A equations he proposed ]2]| . Together 
with our numerical results |2!| , the leading expansion terms obtained in the TB A approach 



read (recall that y = | here) 

e (r) = -— + 0.016781684(2) r 4 / 3 + — ^r 2 lnr - 0.00926790(5) r 2 + C(r 8 / 3 ) . (A.7) 
15 \2'K A 

Now comparing a,2 here with the CPT prediction ( |A.1| ) we obtained k§ = 0.130234474(2) 
of Table 1 (actually this value is deduced from our slightly more accurate results for a<i 
in A 5 _) ). Applying our recipe (|A.4|) to ( |Q|) (in our case = C(i >3 )(i,3)(i,3) = ^ 



|3^1), we then obtain for the coefficient of the r 2 lnr term 

7 3 (^)«^- (A.N) 



12tt 2 27 7 Ks ' 
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Using the numerical value for K5 this relation gives a = 1.49999998(6), which according 
to the "principle of nice numbers" is nothing but |. The numerical value for (3, which can 
be obtained using the as that we found, has not yet been illuminating to us. 

Assuming the validity of ( |A.4j ) with a = | , we can invert the argument and obtain 
a new result. Namely the exact expression for k 5 that we presented (just alluding to the 
"derivation" given here) in 



ls^Sff) 

k 5 = — ^ = 0.13023447336... . (A.9) 

2n 

The same exact expression has been also given recently in [23[ .111 There the authors study 
certain integrable perturbations of the Z^-parafermion CFTs, the case N = 3 correspond- 
ing to where k is same as in A$ . Though conceptually different, their analysis also 
resorts in the last step to some (ad-hoc) prescription for obtaining a finite number out of 
a logarithmically divergent integral — see eq. (32) there. 



12 The K5 read off from eq. (33) of [23] differs from ours by a factor of \/2 due to a different 
normalization of the perturbing field. 
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